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Abstract 

The complete QCD corrections to the scalar top quark pair 
production via 77 fusion in parent e + e~ collider are calcu- 
lated in the Supersymmetric Model. The various origins of 
corrections are discussed. We find that the correction due 
to gluino exchange is smaller than that due to conventional 
gluonic contribution, and the contribution due to squark 
quartic interactions is much smaller. Moreover, the depen- 
dences of the corrections on the masses and the mixing angle 
of the supersymmetric particles are also investigated. 
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1. Introduction 



The Standard Model (SM) |2j has achieved great successes in describing all the avail- 
able experiment data pertaining to the strong, weak and electromagnetic interaction phe- 
nomena up to a few percentage level quantitatively. Only the elementary Higgs boson, 
required strictly by the Standard Model for spontaneous symmetry breaking, remains to be 
found. At present the Supersymmetric extended Model (MSSM)[[J is regarded as one of the 
most attractive model. Apart from containing a large number of experimental consequences, 
the supersymmetric theory is able to solve various theoretical problems, such as the fact 
that SUSY maybe provides an elegant way to construct the huge hierarchy between the 
electroweak symmetry-breaking scale and the grand unification scales. 

One of the most important consequence of the SUSY model is to predict the existence of 
scalar partners to all existing fermions in the SM. In the model two coloured scalar quarks 
(squarks) qi and q~R are required as partners of the standard Dirac quarks. To probe the 
existence of these squarks, to look for pair production qq would be a very promising. The 
pair production may have a comparatively large cross section, depending on the masses of 
the squark. Furthermore, if both pair-produced squarks decay via q — > qj, the signature 
with hard acoplanar quark jets, may be quite good for detection, allowing identification of 
individual events which would not occur in the Standard Model. 

Since q~i and q~R will turn to their mass eigenstates by mixing with each other if they are 
produced, and the mixing size is proportional to the mass of the related ordinary quark H. 
As the top quark is very heavy, the scalar partners of top quark (stop quark) are expected 
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to be strongly mixed. Therefore, one of the two mass eigenstates of stop quarks t% and t 2) 
which are mixtures of the left and right weak eigenstates of the stops q~L and q^, may be 
comparative light (the lightest one of the all squarks) i.e. we may suppose the one of the 
stop mass eigenstates t\ is much lighter than the other one £2- As a result, t\ is very possible 
to be discovered in a relatively low energy range. 

The Next-generation Linear Collider (NLC) is designed to look for the evidences of Higgs 
and other new particles beyond the Standard Model. To detect the existence of the heavy 
squarks, e + e _ collision has advantage over pp or pp collision with its cleaner background. 
Squark pair production produced by e + e~ annihilation has been studied thoroughly, both 
at tree level || and at next-to- leading |J. The overall consideration of producing a scalar 
top pair via 77 collision is still absent and its production rate may be greater than that 
by the e + e~ annihilation as the later has an 's-channel suppression' due to the virtual 
photon propagator when the scalar top quark is so heavy. Moreover, in the process of 
squark pair production via photon fusion, the final particles should not be two different 
squark mass eigenstates, because it is merely electromagnetic interaction phenomenon and 
no weak interaction is involved at the tree-level, other than in the production of e + e~ — > 
QiQj (hJ —1)2) there must be some contribution from Z boson exchange. Therefore, 77 
collision provides another interesting mode of producing squark pair which is worthwhile to 
investigate. Furthermore, with the advent of new collider technique, the NLC is not only 
able to work as e + e~ collider, but also to turn the high energy electron-positron beams into 
the Compton backscattering energetic photon beams with high efficiency J7| . The luminosity 



3 



and energy of the photon collider can be expected to be comparable to those of the primary 
e + e~ collision. Therefore, with the solid support of new experimental techniques, there is 
possibility to yield a scaler top quark pair production directly via the high energy photon 
collision. 

From the above argument we can see that the process of squark pair production via 
photon-photon collisions would be interesting and will be accessible experimentally. We 
know that the study for this process only at the tree-level is far from enough in high pre- 
cision. In this work, we will investigate the 1-loop QCD corrections up to the order 0(a s ), 
including the conventional ones (virtual and real gluon exchange) and those of supersym- 
metric partners. The paper is organized as follows: In Sec. II the analytical formula of tree 
level is presented. The different origins of supersymmetry QCD corrections are discussed 
in Sec. III. In Sec. IV, the numerical results are depicted, and some discussions are given. 
Finally, the conclusions are made. In the Appendix, the form factors appearing in the cross 
section of Sec. Ill are listed. 

II. Tree level result 

In this work, we denote the reaction of stop quark pair production via photon-photon collision 

as 

7(P3, A*)7(P4, V) ► ti(Pl)ti(P2) (2.1) 

where pi and P2 represent the momenta of the outgoing scalar top-quark and its anti-particle, 
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P3 and P4 describe the momenta of the two incoming photons. 

The Feynman diagrams for the process at the tree level are shown in Fig.l-a, and the 
relevant Feynman rules can be found in ||. The corresponding Lorentz invariant matrix 
element of the lowest order for the reaction 77 — > ttfi is written as 

M = Ml + M% + Ml (2.2) 

where Mq, Mq and Mq represent the amplitudes of the t-channel, u-channel and quartic 
coupling diagrams respectively. The explicit expressions can be given as 

4ie 2 2 - 
t — mi 

u — m~ 

Ml = 2ie 2 Q 2 { g^e„(p 3 )e v (p A ) (2.3) 

The Mandelstam variables i and u are defined as i = (pi — P3) 2 , u — (pi—pi) 2 - rrip (i = 1, 2) 
and Qi denote the masses and the charge of scalar top quarks respectively. 
The cross section at tree level can be depicted as 

- r — 1 4vra 2 4 Ar „ ri 8u 4 2u 2 (l - 2u 2 ) , _ , , 

o B \ll - UU] = — Of N C P{1 + + p M (2-4) 

Here fi 2 = mf js, the velocity of the scalar top quarks is defined as (3 = Jl — 4m~ Js. The 
kinematical variable w is defined as w = (1 — /?)/(! + (3). Nq = 3 is the number of colors. 



III. SUSY QCD corrections 
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Since the physical squarks are the mixture of the left- and right-handed weak eigenstates, 
the stop quark mass eigenstates t\ and t% can be expressed by 

t\ = cos 9 1 ti — sm9i tn, i 2 = sin^- tt + cos% Ir (3.1) 

where 9% is the scalar top mixing angle which transforms the scalar top mass eigenstates 
ti (i=l,2) to the weak eigenstates II and £r. We assume that the scalar top quark mass 
eigenstates split heavily, namely mi x « mi 2 . Therefore, in this work we constrain our 
attention only on the calculation of the tit\ pair production. 

The Feynman diagrams for the process 77 — > tit\ including all the QCD corrections up 
to the order 0(a s ) are depicted in Fig.l. The total cross section with the supersymmetric 
QCD corrections corresponding to Fig.l can be denoted as 

a = a B + a soft + 5a {1 ~ loop) (3.2) 

where a so jt and £o-( 1 ~ Zo °p) represent the cross section of soft gluon emission and virtual one- 
loop QCD corrections respectively. 

The one-loop supersymmetric QCD corrections originate from three different sources: the 
conventional virtual gluonic corrections (see Fig.l-b), those due to the gluino exchange and 
those due to squark quartic interactions (see Fig.l-d). In order to show the virtual effects 
from the different sources, we divide the loop corrections into three groups correspondingly. 
In the calculation, the dimensional regularization and on-mass-shell scheme (OMS) || are 
adopted in doing renormalization, whereas when considering the effects from the soft gluons 
infrared divergences being regularized by introducing a small no-zero gluon mass is adopted. 
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3.1 Related self-energy and counterterms. 

The unrenormalized self-energy of the scalar top t\ can be written as the summation of 
three parts. 

E(p 2 ) = SS(P 2 ) + ^~ 9 \ P 2 ) + E«V) (3.3) 

where E^, and E® denote the scalar top self-energy parts contributed by the di- 
agrams the virtual gluon(Fig.l-bl), gluino exchange (Fig. 1-dl) and the scalar top quartic 
interactions(Fig.l-d5) respectively. The squark quartic interactions are involved by intro- 
ducing superpotential of the scalar matter fields [|J. Imposing the OMS renormalization 
conditions ||, one can get the related renormalization constants to 0(a s ) order with the 
formulas 

5M 2 = 6M 2{9) + 5M m + 5M 2 ^ = i?eE(m|), 

5Z = 5Z {9) + 5Z { ' 9) + 8Z® = -i?e^^L 2=m 2 , (3.4) 

dp 1 h 

where Re only takes the real part of the loop integral functions appearing in the self-energies. 
All the definitions of loop integral functions A, B, C and D used in the paper can be found 
in Ref.0. 

The three contribution parts of the scalar top t\ self-energy can be written explicitly as 

zEg(p 2 ) = -i^f {AoKJ + ip 2 {B + B 1 )[p 2 , m g , m h ] + m 2 g B [p 2 , m g , m f j} , (3.5) 

«E (9) (p 2 ) = i-f—d{A [m t } - (mj + m t m~ g sin 20 t ~)-B o [p 2 , m- g , m t ] - p 2 B 1 [p 2 ,m- g ,m t ]}, (3.6) 
iE® (p 3 ) = -ijL{A [m{] cos 2 26~ t + A [m|] sin 2 2%}, (3.7) 
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where m g is denoted as the gluon mass and d = 4 — e is the space-time dimension. 

Note that the self-energy part from stop quartic interaction diagram(Fig.l-d5) expressed 
in Eq.(3.7), is independent on its momentum squared p 2 . This means the corresponding 
part of the squark wave function renormalization constant 6Z® = 0. The other related 
renormalization constants SZ^\ 6Z®, 5M 2( 3\5M 2 ^ and 6M 2 ® can be obtained by using 
the expressions of -Eg. (3.3 ~ 3.7). 

3.2 Renormalized one-loop corrections. 

The renormalized one-loop matrix element involves the contributions from all the SUSY 
QCD one-loop self-energy, vertex, box and quartic interaction diagrams(shown in Fig.l — 
(61 ~ 612) and Fig.l — (dl ~ d7)) and their relevant counterterms. It can be written in the 
form as 

SM (i-ioop) = sm^ + smw + sm® 

(3.8) 

= ie 2 Ql eM^M [fi <T + h P?Pi + fs p1p u 2 + h P2P1 + h T%&\, 

with 

h = f% + f'f ) + ff ) , (3.9) 

where SM^f} represents the gluon correction part of the one-loop renormalized amplitude, 
and 5M^\ SM® represent the gluino exchange and the squark quartic interactions con- 
tribution respectively. All the renormalized self-energy, vertex, box and quartic interaction 
loop corrections due to the different sources, are summarized in Lorentz invariant factors 
fvil,ii fi^ an d fi^ (i = 1, 5), respectively. By using this formulation we are to show a clear 

8 



distinction among those three different origins of supersymmetry QCD corrections. The 
explicit expressions of form factors can be found in Appendix. 

Now we can achieve the one-loop SUSY QCD corrections of the cross section for this 
subprocess in unpolarized photon collisions. 



5a {1 - loop \s) = — dt 2Re E (A/J • 5M^ ioop ») 

1 * spins 

(3.10) 

where = (m| — |s) ± \s(3. The bar over the sum recalls averaging over initial spins. The 
correction 5<t^+*) = + will be used to denote the corrections arising from the pure 
supersymmetric particle contribution. 



3.3 Real gluon emission 

The real gluon emission diagrams for the process 77 — > titig are drawn in Fig.l — cl ~ c7. 
The cross section of real gluon emission is deduced as 

N c 1 f d 3 p l( l 3 p 2 d 3 k 4 _ 



with 



E l^e a /i 2 = e (fL ■ + o • |M«| 2 + /*, • 

spins spins 

+ 2ft al ■ Rc(MqMq) + 2fti ■ Re(M^Mi) + 2f% al • Re(M^M^)) 



(3.116) 

— * 

where k and to denote the momentum and the energy of the outgoing gluon respectively, and 
Mq (a = t, u, q) are the invariant matrix elements of tree level defined above. The factors 
fre a i (a, 6 = t,u, q) can be found in Appendix. 

In the region around k ~ 0, the cross section in Eq.(3.11), corresponding to the cross 
section for a soft gluon being emitted, may be approximately alternated as 

where 

h = 2vr hlog[—\ + j (3.13a) 

1 - 2u 2 f , . ,2A£, . 23 . 1, 9r .] 

hj = 2vr— ^- hlog[w]log[—} + 2^2!^] + -log 2 [w]\ (3.136) 

Here a small soft gluon energy cut AE is introduced. Li2\x\ is the Spence function. The 
definition of 3 is the same as in Eq.(2.4). The expressions Eq.(3.12) and Eq.(3.13a,b) are 
similar to those as described in Ref.|| for a soft photon being emitted, that just is a check 
of Eq.(3.11) for the real gluon emission in a certain sense. However, for the cancellation of 
the infrared divergences due to gluons, it is sufficient that only the soft gluons , i.e. gluons 
with energy uj < AE, are relevant, where AE is small, which may be set to be smaller 
than all relevant energy scales in concerned problem. Therefore, we adopt the soft gluon 
approximation method with a cut AE in numerical calculations for the soft gluon emission 
corrections. With this approach the radiation from internal colour lines or quartic vertices 
does not lead to IR-singularities and can be neglected in this approximation. 
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The correction to the total cross section from the gluonic QCD contributions should be 
expressed as below 



6&te) = 5a%l + a aoft (3.14) 

IV. Numerical results and discussions 

In the numerical evaluation, we take the value of the top-quark mass to be m t = 175 GeV, 
a = 1/128 and use the one-loop running coupling constant a s . We assume < mi 2 for 



the masses of the stop quark mass eigenstates. The present constraint from experiments flO 
is that the lighter scalar top mass eigenstate mi x should be heavier than 50 GeV. 

The corrections to the cross section of the process 77 — » t\t\ as a function of the center- 
of-mass \/I are studied with the values of the involved SUSY parameters given below 

m k = 200 GeV, m h = 550 GeV, m- g = 250 GeV, cos0~ t = 0.7 (4.1) 

The results are depicted in Fig. 2 (a). These results of numerical calculations show (i) Both 
corrections 5a^ and 5a^ +t ^ reach their maximal values around the threshold of titi pair 
production and vary from positive to negative value as increasing Vs. The former varies 
severely and in general is several times higher than the later, (ii) The relative corrections 
due to the gluon contribution can even reach +50% when c.m.s energy \/~§ is just above the 
pair threshold, whereas the maximal value of Sa^ +t ^ /&b is only about +10%. Therefore, one 
may see the fact that the SUSY QCD corrections are dominated by the gluon contributions. 
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The t\ti pair production via photon-photon fusion is only a subprocess of the parent 
e + e~ linear collider. It is easy to obtain the total cross section for the scalar top quark pair 
production via photon fusion in the e + e~ collider, by simply folding the cross section of the 
subprocess (7(77 — > Mi) with the photon luminosity. 

a(s) = / dz— —a(jj -»• at s = z 2 s), (4.2) 
J 2™,^ / yfs az 

where yfs and vl are the e + e~ and 77 center-of-mass energies respectively and dL^/dz is 
the distribution of photon luminosity, which is expressed as 

^ = 2 2 P" (4.3) 

LIZ J Z /Xmax 

where / 7 / e is the photon structure function of the electron beam |7|)[lT||. We take the structure 



function of the photon produced by the most promising Compton backscattering as [0, |12 



fComp= \ ^(l-^ + i^-^y + s^), forx < 0.83,^ = 2(1 + ^) 
7/6 \ 0, forx > 0.83. 

(4.4) 

The cross sections of the process e + e~ — > 77 — > t±ti versus \/s, with the same condition of 
Eq.(4.1), are depicted in Fig.2(b). There the magnitude of the corrections is also significant 
for this process and the sign of corrections goes from positive to negative with the increasing 
of y/s. The total cross section involving SUSY QCD corrections shows the weak dependence 
on the parent c.m.s energy when ^/s > 1 TeV. Furthermore, by analysing our numerical 
data, the whole correction 5a/ as can reach +35%. 

In Fig. 3 the concerned corrections are plotted as functions of the light scalar top mass 
m^p with = 1000 GeV, m^ 2 = 550 GeV, rrig = 250 GeV and cosOi = 0.7. In most regions 
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of the value m^, the corrections due to gluonic contributions are much larger than those 
due to the scalar top quarks and gluino contributions. Around the region of ~ 500 GeV 
where just reaches its pair production threshold, the gluonic corrections become positive, 
whereas they are negative when y/I is far above the threshold. 

In Fig. 4 the concerned corrections with various m t ~ 2 and = 500 GeV, m~ tx = 200 GeV, 
nig = 250 GeV, cosQi = 0.7 as well, are depicted. Since the heavy stop mass only appears 
in the amplitudes of the diagrams involving the squark quartic vertex, the dependence of the 
corrections on mf 2 manifests only in the correction Sa^\ The correction 8a^ /a B approaches 
to zero with the increment of mi 2 and the maximal correction is only about —1%. Note that 
when m t ~ 2 = 550 GeV , the relative correction 5a^/a B is only —0.2%, whereas in Fig. 2(a) 
the correction 8a^ t+ ^ / a B with the same condition is about +4%. That shows 6 a® is always 
smaller than both Sa^ and 8a^\ Moreover, the curve of the correction due to the squark 
quartic interactions shown in Fig. 4, becomes rather flat in the large mi 2 region, which implies 
the dependence of 6 a® on heavy t 2 is very weak. 

The gluino exchange contribution to the cross section is calculated as a function of 
nig. The results can be seen in Fig. 5, where the parameters are set as = 500 GeV, 
= 200 GeV, m~ t2 = 550 GeV and cosOi = 0.7. The curve of the correction 5a^/a B is 
nearly parallel to that of the gluon and stop contribution 5a^ 9+£) /a B , which reveals that the 
correction dedicated by the gluino exchange contribution is independent on m g , when the 
mass of gluino becomes heavy. 

The dependence of the corrections 5a^ +t ^ on the mixing angle cosOi f° r = 500 GeV , 
= 200 GeV , mi 2 = 550 GeV and m g = 250 GeV is shown in Fig. 6. The combined cor- 

13 



rection Sa^ +t ^/crB varies from —3% to +5%. From the figure one can see that the correction 
will reach its peak values when cosOi = ±0.7, namely when the SUSY weak eigenstates ti 
and i R are mixed equivalently. 

VI. Conclusion 

In this work we have calculated the complete 0(a s ) order SUSY QCD corrections, includ- 
ing the conventional gluon contribution and the pure SUSY particle contributions as well, 
to the tjti pair production process via photon-photon fusion in the initial e + e~ collider. 

From the results of numerical evaluation, it is expected that the corrections due to the 
pure SUSY particle U {i = 1, 2) and g contributions are smaller than those from the virtual 
gluon loop diagrams. Among the pure SUSY particle contributions, those due to the squark 
quartic interactions are much smaller than those due to the gluino exchange with reasonable 
parameters. The gluonic correction may reach +50% in the photon fusion subprocess, when 
VH is around the t\ti pair production threshold, whereas the corrections there contributed 
by ti and g diagrams are only about +10%. The total QCD corrections to the subprocess 
will be most significant when is around the pair production threshold. For the cross 
section of the process e + e~ — > 77 — > t±ti, the SUSY QCD correction at one- loop level can 
reach 35%. Furthermore, the corrections also show the decoupling with the heavy g and £2 
masses, and the dependence of the QCD corrections on the mixing angle d~ t obviously. 

One of the authors, Z.H. Yu, would like to thank the Institute of Theoretical Physics 
of University Vienna for the warm hospitality extended to him during his stay under the 
agreement the exchange program (project number: IV.B.12). 
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Appendix 



The form factors f- 9 ^ due to the gluon contribution can be expressed as 

/i>,i = S ^f{ B o\pi+P2,m il ,m il }- 

((m 2 g C + m\2C xl ) - 2(p! • p 2 )(2C + C n )) 
[Pi, -Pi -P2,m g , m~ h , m t -J } 

- ^r 1 {(#o[Pi - p 3 , "V ™tj - 2C 24 [pi, -P3, m 9 , w^, mfj + CWf-Pa, ~P4, m^, m^, m t ~J + 
(2(p! • p 2 )(2D 27 + ^n) + 2( Pl - p 2 ) • p 3 (Au2 - Am) - ^27 - 2m? £311) 
[Pi, -P3, -P4, m s , m fl , m fi , m t -J) + (p 3 <-> p 4 )} 

+ 28Z<*) 

(Al) 



/ife = - £ &{(2C , + 3C 11 -C 12 + C 21 -C 23 )[pi ! -P3,m s ,m fl ,m ?1 ] 
-(C12 + C 23 )[-p 3 , -p 4 , m h ,m h ,m h ] 
+ (2(D 312 + D 313 - D 27 - 2D 311 ) + 

m 2 g (D u - D 12 + D 21 - D 24 - D 25 + L> 26 ) + 

2m| (D 21 - D 24 + A31 + A310 - ^34 - ^35) + 

2(pi • P2) (^34 + D35 - 2£>ii + 2D 12 - 3D 21 + 3^24+ (A2) 
2D 25 - 2D 26 - D 31 - D 310 ) + 

2(Pl - P2) ' P3(£>22 - £>24 + £>25 ~ £>26 ~ £>34 + D 35 + 

Ahj - D 37 - D 38 + D 39 )) 
[Pi , -P3 , -Pa , m g , m h , m fl , m h } } 
+ (Pa «-»■ P4) 
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/& = 9 ^ Tr ^{A [m il ] + (mlB + 4t-(B + B 1 ))[p 1 -p 3 ,m g ,m il ]} 

2 

l^t^f {( 2B ° + m g , m h \ + (2B + Si)[pi - p 3 , m g , m h ] + 

H (m 2 g (C + C n ) + (C + 2C n + C 21 ) - 4C 24 - 
2(pi • p 3 )(2C + 3C n + 2C 12 + C 21 + 2C 23 )) 
\pi,-P3,m g ,m~ tl ,m il \^ 

+ -V<5^) 

t—mt 
H 

+ {2(C 12 + C 23 )[pi, -p 3 , m 9 , m fl , m t -J + 

2(Ci2 + C 23 )[-p 3 , -p 4 , m h ,m k ,m k \- 
(2(2D 27 + L> 311 + .D,!, - L> 313 )+ 
m 2 g (D 26 - D - D u - D 12 + D 13 - D 24 ) + 

2( Pl ■ p 2 )(2D + 3D n + 2D 12 - 2D 13 + D 21 + 3D 24 - 

D 25 -2D 26 -D 310 + D U ) + 
2 (pi - p-i) • Pz{D 12 - D 13 + D 22 + D 23 + D 2A - D 25 - 

2D 26 - D 310 + D 36 - D 38 + D 3g )) 

[Pi , -Ps , -Pa , m g , m h , m h > m *i 1 ~ 
(2(D 3 n - + D 3l3 ) + m2(L> 26 - D 25 ) + 2m\{D 3W - D 35 ) + 

2(pi • pa) (2^25 - 2D 26 - ^310 + ^35) + 
2(P1 - P2) • P4(-D310 - £> 37 - £>38 + £39)) 

[pi , -p 4 , -P3 , m g , m fl , m h , m h } } 

fit A = fi 9 \i^u, V3 ^ Vi ) 



/& = - £ &{(2C + 3C 11 -C 1 2 + C2i-C23)[pi,-P3,m 9 ,m ?1 ,m ?1 ] 
-(C12 - C 23 ) [-p 3 , -p 4 , m h , m h , m fl ] 

+(m 2 g (D 13 + L> 26 ) + 2771^(^25 + D 310 ) - 2(D 27 + £> 312 + L> 3 i 3 ) + 

2(Pl - P2) • P3(0>23 - 0> 26 - £>38 + D39)- 

2( Pl -p 2 )(2D 13 + D 25 + 2D 26 + D 310 )) 
[Pi , ~P3 , -Pa , m s , m fl , m h , m fl ] } 
+ (P3 <-> P4) 
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The form factors f- due to the gluino exchange can be expressed as 

- 2{d • C 24 + 2m?C + 2m}C l2 + 2( Pl ■ p 2 )(C + C 12 ) + 

2 (P2 • P3)(C + Cn) + 2(pi • p 3 )Ci2 + m t m- g s%n2Q~ t C ) 
[-P3, -P4,m t ,m t ,m t ] 

- 4(m 2 g (m 2 - m?)A) - 2(m 2 + m?)D 27 + 

m\{m 2 D X3 + mf(£> 13 - A) - 2Ai) + m| i (2£) 23 - D 13 - 2D 25 ))+ 
(Pi ■ P2)(m 2 t (D 13 - D u ) - AD 27 + mf (Ai + D 13 ) + 

m?(Ai - D 13 + 2D 2l + 4D 23 - 4D 25 )) + 
(Pi • P 3 )K(Ai - As) + ™1(A + Ai + 2D 12 - 3D 13 ) + 

m 2 h {D 13 - 2D 23 + 2D 26 )) + 
(p 2 • p 3 )(m 2 (D 12 - 2D l3 ) - mf (A + D 12 ) + 

mf (2Di 3 - - A 2 - 2£ 23 - 2D 24 + 2D 25 + 2Ae))+ 
2( Pl ■ p 2 y(D 23 - As) + 2(pi • p 2 )(p 2 • p 3 )(A 6 - As) + 

2 (P1 -P2)(P1 -P3)(D 25 + Ae - A3 - Am) + 

2(pi • Pa)(P2 • Ps)(D 12 - D 13 + D 22 + D 23 - 2D 26 )+ 
m t m~ g sin26~ t ((m 2 t - m 2 g )D - AD 21 + m\(2D 13 -D - 2D n )+ 
2(pi • p 2 ) A3 - (P2 • Ps)(A + 2^13) + 
(Pi • Pa) (A + 2Ai + 2D 12 - 4L> 13 ))) 
[Pi, ~P3, -P4, m 5 , mt, m t , m t ]} 
+ (p 3 ^P4) + 2 5Z«) 



= 9 -£f{(Ci 2 + 2C 23 )[-p 3 ,-p,,m t ,m t ,m t ] 

- (m 2 (Ai - 2D 12 + D 13 + 2D 21 - 2D 24 - 2D 25 + 2D 26 ) + 

m?(A + 3£>n - 2D 12 - D 13 + 2D 21 - 2D 24 - 2D 25 + 2D 26 ) + 
ml(D 13 -2D 23 + 2D 25 )+ 

2((pi • p 2 ) - (P2 ■ Pa)) (As - D 12 - D 23 - D 24 + As + ,026) + 
m t m g sin29i{D Q + 4(Ai - A2 + Ai - A 4 - As + Ae))) 
[Pi> -P3, -p4,TOp,m t ,m t ,m t ]} 
+ (P3 <-+ P4) 
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f® = J'i-% )3 {(K + m t m 9 sin29i)B Q + t ■ Bjfa - p 3 , m- g , m t \ - A [m t }} 
+ Jf- m \ ) { rf ' 5 o[-Ps, mt, m*] + 4 • (2( Pl ,p 3 )C 12 - m?Cn- 

m? t -(2C + C 11 )-m|C 11 - 
2m t m § sm20 t ~ (C + Cn))[pi, -p 3 , rn,g, m t , m t }} 
^(t ■ 6Z® - SM^ - ml ■ 8Z<*>) + j+rSZ® 



+ ^5§r {(C + C n + 3C 12 + 4C 23 )[-p 3 , ~P4, m t , m t , m t \ 

- (2m 2 t (D 13 - D n - D 12 - D 24 + D 26 ) - m|L> n + 
mf (2(L> 13 - - L>i 2 - D 24 + L> 26 ) - 3D ) + 

2{p 1 -p 2 )(D 13 -D 12 -D 24 + D 26 )+2{p 1 -p 3 )(D 12 -D 13 ) + 1 8j 

m t m 5 sin29i (4(L> 13 - D u - D 12 - D 24 + D 26 ) - 3D )) 

\pi, -P3, -P4, rrig, mt, m, m] 

- (mj(Do - D n + L> 12 - D 13 - 2D 25 + 2^20) + 4D 27 + 
m 2 t (D u - D 12 + D 13 - 2D 25 + 2D 26 ) + 

ml (D l2 - D 13 - 2D 21 - 2D 23 + 2D 24 + 2D 25 - 2D 26 ) + 
2(pi • P2)(D 25 - D 23 ) + 2{p 2 ■ p 4 )(D 23 - D 26 )+ 
2( Pl ■ p 4 )(D 13 - D 12 - D 22 + D 24 - D 25 + D 26 ) + 
m t m~ g sin26~ t (D - 4L> 25 + 4L> 26 )) 
[pi, -p 4 , -p 3 ,m~ g ,m t ,m t ,m t ]} 

ff ) = ff ) (i^u,p 3 ~p 4 ) (A9) 

ff = 9 -^f{(C + C 11 + 2C l2 + 2C 23 )[-p 3 ,-p 4 ,m t ,m t ,m t } 
- (ml(D u + D 12 - 2D 13 + 2D 24 - 2D 26 ) + 
m 2 ,(D + D 12 + 2D 13 + 2D 26 )+ 

m 2 (2D 13 - D 12 + 2D 26 )+ (MQ) 
2(pi • p 3 )(D 13 - D 12 - D 22 + D 26 ) + 
m t m~ g sin26~ t (D Q + AD 13 + 4L> 26 )) 
[Pu -P3, -PA,m~g,m t ,m t ,m t \} 
+ (P3 <-»• Pi) 

The factors arising from the contribution of squark quartic interactions, are given as 
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follow 



f? = -^2Cos 2 26i {4 • C 2 i[-p3,P3 + P4, mi, m in m h\ ~ B o[Pi + P2, ra^, m t -j} 
+ (cos — > sin, — > m t ~ 2 ) 



/ 2 = fftcos 2 2^- (C 23 - C 22 ) [-p 3 , Ps + P4, m h , m fl , m t -J , Al2 \ 
+ (cos — > sin, — > m^) 



/ 3 W = ff!cos 2 2# t ~ (C 23 - C22) [-P3, Ps + Pa, m h , m h , m t ~J 
+ (cos — > sin, — > m^) 



/i i) = /i*" ) (*-fi) (^14) 



if = / 2 W (^15) 
The factors shown in Eq.(3.11b) have the expressions as 

ffeai = 2./ 2 -mf 1 .(/ 11 + / 22 ) + 2(p 1 .p 2 )./ 12 

(Pi -P2) 1 2 (Pi - Pa) 1 (Pi - Ps) 2 



(Pi • P2) 2 (Pi • P3) 2 (Pi • P2XP1 • Ps) (Pi • P2XP1 • P3) 



(AW) 



fZl = frUP3^P*) (A17) 
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fZ l = 2(p 1 -p 2 )I 12 -mUl 11 + L 



122 



(A18) 



2 • I 2 - m? ■ (I u + I22) + 2(pi • p 2 ) ■ I12 

2m? 
ti 

(Pi • P2) 
2m? 



9 

(/! + J 2 ) + -;—^h + 



7?) 



(Pi-Ps) (Pi-Pa) 
2 



m h T , (P3 - Pi) • P2 r . (P4 - Pi) • P2 r 
Ji H ; ^ ±2 1 ; ; 



(Pi ■ Ps) 



(Pi • P4) 



m 2 



(pi • P2) 2 (Pi • P2) (Pi • Ps) (Pi • P2) (Pi • Pa) 

(P2 • Ps) (P2 • Pa) 



ra- 
ti 



+ 



+ 



(Ps • Pa) 



(Pi • Ps) (Pi • Pa) (Pi • P2) (Pi • Ps) (Pi • P2) (Pi • Pa) (pi • Ps) (Pi • Pa) 



ff eal = I 2 - m 2 r ■ (In + I 22 ) + 2(p! • p 2 ) ■ h 2 



+ 



m~ 



-h - 



2 

m| 
ti 



h- 



m~ 



(Pi • Ps) (Pi • P2) (Pi • P2) 



T , (Ps - Pi) • P2 T 

I2 H 7 v — H 

(Pi • Ps) 



0419) 



(420) 



freal ~ f?Ll(P3 ~> Pa) 



(A21) 



lab.. 



(Pa ■ k)(p b ■ k)... 



(a, 6 = 1,2) 



(A22) 
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Figure Captions 

Fig.l Feynman diagrams including the contributions from the tree-level and next-to- 
leading order in Supersymmetric QCD for 77 — > t±ti process: (1-a) tree level diagrams; 
(1-b) corrections due to gluon contribution; (1-c) soft gluon emission; (1-d) corrections due 
to the pure supersymmetric particle g and U. The diagrams by exchanging initial photons 
are not shown here. 

Fig. 2 When m~ g = 250 GeV, m k = 200 GeV, m h = 550 GeV and cos6~ t = 0.7: (a) 
the cross section 17(77 ~~ Mi) as a function of V§, solid-line for tree-level contribution erg, 
dotted-line for the conventional QCD correction a^ 9 \ dashed-line for the pure supersym- 
metric particle QCD contribution cr^+*); (b) the cross section a(e + e~ — * 77 — > tit±) as a 
function of y/s, solid-line for tree level contribution a B , dash-dotted-line for the overall QCD 
correction a^ 9+9+t \ 
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Fig. 3 The cross section a as a function of m^, when = 1000 GeV , m ( - 2 = 550 GeV ' , 
mj = 250 GeK and cosQi = 0.7: solid-line for tree-level contribution a B , dotted-line for the 
conventional QCD correction o - ^, dashed-line for the pure supersymmetric particle QCD 
contribution a^ +t \ 

Fig. 4 The relative correction 8a/a B as a function of m t ~ 2 , when = 500 GeV, = 
200 GeV, m- g = 250 GeV and cos0 f = 0.7: dotted-line for the contribution 8a^ 9+ ^/a B , 
dashed-line for the QCD correction due to squark quartic interactions 8a® /o B . 

Fig. 5 The relative correction 8a ja B as a function of m g , when V§ = 500 GeV , = 
200 GeV, mi 2 = 550 GeV and cosOi = 0.7: dotted-line for the contribution 8a^ 9+t ^/a B , 
dashed-line for the QCD correction due to the gluino exchange 8a^/a B . 

Fig. 6 The relative correction 8a /a B as a function of cosOi, when y/J = 500 GeV, = 
200 GeV, mi 2 = 550 GeV and m g = 250 GeV: dotted-line for the gluon contribution 
8a^/a B , dash-line for the pure supersymmetric particle contribution 8a^ +t ^ /a B . 
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Fig.l (continued) 



vAAA^r \AAA-j q 

% 'I 

I % |f 

\A/VO— - vAAA^ 

l-b9 1-blO 



I I I ° 
vAA/ 1 '- vAAAJ 

l-bl2 1-cl 



vAAAj^ vAAAp 

i%. i 

vAAA> vAAAJ- 

l-c3 l-c4 



Fig.l (continued) 



Fig.2(b) 
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